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ABSTRACT 

We investigate the role of neutron star superfluidity for magnetar oscillations. Using a 
plane- wave analysis we estimate the effects of a neutron superfluid in the elastic crust 
region. We demonstrate that the superfluid imprint is likely to be more significant 
than the effects of the crustal magnetic field. We also consider the region immediately 
beneath the crust, where superfluid neutrons are thought to coexist with a type II 
proton superconductor. Since the magnetic field in the latter is carried by an array of 
fluxtubes, the dynamics of this region differs from standard magnetohydrodynamics. 
We show that the presence of the neutron superfluid (again) leaves a clear imprint on 
the oscillations of the system. Taken together, our estimates show that the superfluid 
components cannot be ignored in efforts to carry out "magnetar seismology". This 
increases the level of complexity of the modelling problem, but also points to the 
possibility of using observations to probe the superfluid nature of supranuclear matter. 



1 INTRODUCTION 



Anomalous X-ray pulsars (AXPs) and Soft Gamma-ray 
Repeaters (SGRs) are widely believed to be magnetars; 
neu tron stars powered by an ultrastrong magnetic field 
(see I Woods fc Thompsoiil (|2004l ) for a review) . Observations 
(mainly conducted by X-ray satellites) have established ba- 
sic parameters like the magnetic field intensity, B ~ 10^^ G, 
and spin period, P ~ 10s, for this class of objects. They 
have also revealed a complex emission pattern with alter- 
nating periods of burst activity and quiescence. SGRs, which 
are typically more active than AXPs, are the only ones ex- 
hibiting giant flares. These flares, which are believed to be 
triggered by some sort of in s tability in the magnetic fiel d 
l|Duncan fc ThompsonI Il992l : I Thomson fc DuncanI Il995l ). 
are by far the most energetic events associated with magne- 
tars. 

An exciting contribution to magnetar phenomenology 
was provided by the recent discovery of quasi-periodic os- 
cillations (Q POs) in the late tail spectrum of the two 
giant flares (llsrael et all l2005l : IStrohmaver fc Wattd l2005l : 
IWatts fc StrohmaveiibOOg T There may also be evidence for 
a single QPO in the d ata of the thir d known flare, observed 
back in the 1970s ( Ba rat et aj|l983 ). The frequencies of the 
most prominent QPOs lie in the range 30 — 100 Hz, exactly 
where one would expect to find the seismic osci llation modes 
of the magnet ar's crust (jHansen fc Ciofi]||l980l ). This is con- 
sistent with the theoretical expectation that the energy re- 
leased in a giant fiare is sufficient to fracture the crust and 
excite its normal modes (lDuncarJl998l ). If this interpretation 
of the QPOs is correct then we may have the opportunity 
to carry out magnetar " asteroseismology" ; a comparison be- 
tween theoretically predicted mode frequencies and the QPO 



data, with the ultimate goal of constraining the prope rties 
of neutron star matter ijSamuelsson fc Anderssonll2007l ). 

Indeed, several recent papers have attempted to con- 
strain the bulk equation of state of neutron star matter, 
assuming that the observed QPO frequencies are identi- 
fied with the fir s t few toroidal seismic modes of the crust 
(see IStrohmaved l|2007l ) for a recent review and references). 
This is natural as a first step, but in reality the sit uation 
is like ly to be more complicated. As suggested by 
(2006") and ^Glam pedakis. Samuelsson fc AnderssonI ( 
the strong magnetic field would likely couple an oscillat- 
ing crust to the liquid core on a very short timescale. Then 
the observed QPOs would be the manifestation of the cou- 
pled crust-core dynamics rather than the dynamics of the 
crust alone. Possible evidence that the magnetic core plays 
an active role is given by the presence of a low frequency 
QPO in the data of the December 2004 flare in SGR 1806- 
20. This QPO is d ifficult to reco ncile with the seismic mode 
interpretation (Isr ael et "ail l2005l ). It is therefore conceivable 
that magnetar "seismology" also probes the (much less well 
known) properties of the interior magnetic field. This obvi- 
ously comes at a price. We now have to model global crust- 
core oscillations, a problem that is considerably more chal- 
lenging than that of pure seismic crust modes. 

Another aspect of neutron star physics, which is directly 
relevant to the mode interpretation of the QPOs, has re- 
ceived almost no attention so far. Young and mature neutron 
stars (older than a month or so) are sufficiently cold that the 
bulk of their interior liquid matter is in a superfiuid state. 
In the crust, for densities above the nuclear drip density 
p ~ 4 X 10^^ g/cm^, the "dripped" neutrons are expected 
to form a superfluid below a temperature ~ 5 x 10^ K. The 
dynamical role of these "free" neutrons could be important. 
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After all, they account for ~ 80% of the total mass in the 
crust. Similarly, in the liquid core we expect to find both 
neutrons and protons in a superfiuid state (below a simi- 
lar threshold temperature). In the outer core, the protons 
most likely form a type II superconductor, provided that 
the interior magnetic field does not exceed a critical value 
« 3 X 10^® G. As a consequence, any magnetic field that pen- 
etrates the proton plasma will form a large number of q uan- 
tised magnetic fluxtubes (jBavm. Pethick fc Pineslll969l ). 

It is clearly relevant to ask to what extent the physi- 
cal components (the crust and the magnetic field) that play 
the leading role in the magnetar QPO problem are sensitive 
to neutron and proton superfluidity/superconductivity. The 
aim of this investigation is to provide an insight into this 
issue, and improve our understanding of the relative impor- 
tance of the multifluid aspects of the magnetar oscillation 
problem. By carrying out a local analysis, i.e. considering 
uniform parameter models, we learn how the shear waves 
in the neutron star crust are affected by the presence of a 
superfiuid neutron component. Similarly, a local analysis in 
the core tells us how the Alfven waves are altered by the 
presence of the neutron superfiuid (which provides the bulk 
of the core mass). Not surprisingly, the entrainment between 
neutrons and protons turns out to be the key parameter in 
these problems. This initial (order of magnitude) analysis 
serves as a useful guideline for future (more detailed) work 
for reahstic neutron star models. 



should be noted that we are not accounting for effects due 
to the presence of neutron vortices, e.g. the mutual friction 
and the vortex tension here. In principle, these effects will 
be present even in the very slowly rotating magnetars, and 
it will be interesting to consider them at a later stage. Our 
initial aim is to explore the leading order effects of this rather 
complicated problem. 

The corresponding equation of motion for the crust nu- 
clei can, labeling the variables associated with the nuclei 
with the constituent index c, be written 

(1 - ec)9?e. + ecdU^ + + V, V,<I> 

- (V.C^)V,/ic + V.AcAc = Ac/r' + Ac/r^ (3) 

Here Ac representing the Lagrangian variation along the 
crust motion (associated with a displacement Q.), and 



Ac/ic 



V,(ncCc 



+ (|^) V.(n„a+CcV,/ic (4) 



It is also worth noting that 



(5) 



The charged component equation includes both elastic 
and magnetic contributions. The former can be written 



(6) 



2 MULTIFLUID DYNAMICS OF THE CRUST 

2.1 Lagrangian perturbation equations 

We want to model linear perturbations in a neutron star 
crust penetrated by a superfiuid neutron component. It 
is natural to use a Lagrangian picture to describe this 
problem. Hence, we c ombine the two-fiuid Lagrangian per - 
turbation equations l|Andersson. Comer fc GrosartI |2004| ) 
with the relevant elastic terms and the magnetic force 
iGlampedakis fc AnderssonI 120071 ). Since all known magne- 
tars are slowly rotating, with periods of several seconds, it 
makes sense to focus on the non-rotating problem. Then we 
need an equation for the superfiuid neutron displacement 
which can be written 

(1 - e„)9?ff + £ndk" + V«cS$ + C^V, V.cE. 

- (V,Cn)Vj/in + V,An^in = (1) 

We label the variables associated with the neutrons with a 
constituent index n. An represents a Lagrangian variation 
along the neutron fiow (associated with the displacement 
5n). The variable Sn (assumed constant here) encodes the 
entrainment effect, 5^ represents the (Eulerian) perturba- 
tion of the gravitational potential "1> and jln is the (specific) 
chemical potential for the neutrons. We have 



where 



An/in 



V,(nnCn) 



+ (|^) V,(ncec)+C;V,Mn (2) 

^ / Tin 

where and ric are the number densities of the neutrons 
and the baryons making up the crust nuclei, respectively. It 



(7) 



(here one should not confuse the shear modulus jj, with the 
chemical potentials ^x) Meanwhile, the magnetic term fol- 
lows from the standard electromagnetic Lorentz force. That 
is, in this case we have 



rmag _ rL _ J_ jj r^k 



(8) 



Eliminating the total current with the help of Ampere's law 
(i.e. r = (c/47r)e'-''°VjBfc), this becomes 



47rpc 



(9) 



Wor king out the L agrangian variation using 
(Glampedakis fc Anderssonl l2007) 



we arrive at 



and 



A, ( — 1=0 



(10) 



(11) 



AcB, = B^V^ii - B,:(V,e^) + V,e> (12) 
Finally, using 

Ac(V,BO = Vj(AcBO - BiV.V.^i (13) 
we obtain from (|9]) 



B' 



■[V,(AcBO- V.(AcB,)] (14) 
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These are all the relations we need to solve the problem. 
As far as we are aware, this is the first time that the per- 
turbation problem for combined superfluidity, elasticity and 
magnetic flelds has been formulated. The equations we have 
given can be directly applied to studies of global mode os- 
cillations of a superfluid neutron star with a crust. 

For later convenience it is useful to note that we could 
equally well have worked with Eulerian variations. In fact, 
since /l = in the background configuration we must have 
Ac/™''^ = 5/™''^. Moreover, one can show that in the case 
of an incompressible fluid and a uniform background field 
(see below) we have 



(15) 



where "hats" denote unit vectors. We have also defined the 
Alfven wave velocity 



2 
"A 



47rpc 



(16) 



It is important to note that in the superfiuid system the 
Alfven velocity scales with the number density of charged 
nucleons, not the total baryon number density l|Mendelll 
Il998f ). 



2.2 Plane- wave analysis 

As a first step towards understanding the problem, let us 
consider the simple case of a uniform, non-rotating back- 
ground. For an incompressible model, we have 



V.^i = — * Ax/ix = 



(17) 



Then the problem simplifies to (note that we will have 5$ ■ 
Vi<I> — for a uniform background) 



(1 ~ en)dUl + £ndU^ = 
(1 - e.)dk^ + £cdU? = Ac/r' + Ac/r 



We also need 



(18) 



(19) 



(20) 



(21) 



where the shear velocity, Vs, is defined by 

2 M 
Vs = — 
Pc 

Note that the shear velocity scales with the number density 
of nucleons locked in the crust lattice, not the total nucleon 
number density as would be the case if there were no su- 
perfluid component. This distinction has not been made in 
previous work where th e crust is mode ll ed as a sing le com- 
ponent, see for example [Puncanl (|l998l ): IPird l|2005l l. 

We now consider (short wavelength <^ the radius of the 
star) wave-propagation in this system. Making the st a ndard 
plane-wave Ansatz (see ISiderv. Andersson fc Comeij (|2008l ) 
for a recent analysis of the analogous non-magnetic two-fluid 
problem) 



where the index x is either n or c, we have 
fc'Af = 



(22) 



(23) 



i.e., the waves are transverse, and 



Ac/f 



-vie [{B^ 



■ k,{Bjec){B'k 



From Eq. (|18|l we immediately get the relation 

AC 



1 



(24) 



(25) 



Using this in Eq. (I19p we arrive at 



v'ik'iB.k'f 



-vlk'{Bjy){&A'i)h (26) 



where we have introduced 



' En 



'En £c 



Deflning 



2 2,2 



(27) 



(28) 



the frequency of "pure" elastic waves, and the Alfven-wave 
frequency 



2 



we have an equation for Al 

,2 

{Bjk-yu}\ 



2 7 2 

v^k 



(29) 



2 



A^, 



~iol{B,y){B' A'i)k (30) 



In order to arrive at the final dispersion relation we first 
note that contracting the above equation with fc' leads to 
the constraint 



Lol{Biy){&A'i)=Q 



(31) 



Thus, we can either choose to look for solutions where fc' 
is orthogonal to the local magnetic field or we see that the 
polarisation Al, and hence A", must be orthogonal to both 
fc' and B^ . Since the right hand side of (|30|l vanishes in both 
cases we find that all non-trivial solutions must be such that 

,2 



tJo 



{Bj^fiul = 



That is, we have the general dispersion relation 

, ,2 ■ 



2 



{B.k^r 



(32) 



(33) 



Note that, in the degenerate case when Bik^ — we cannot 
uniquely determine the polarisation; it can lie in any direc- 
tion in the plane orthogonal to fc*. Also, it is clear that such 
waves do not depend on the magnetic field at all. Generi- 
cally the polarisation is, however, well defined (up to scale 
since we have a homogeneous system) to be orthogonal to 
both k' and B\ 

We need to estimate the magnitude of the various terms. 
Let us first focus on the entrainment. Using ^ we find that 



= 1 - — 

Xc 



(34) 



where Xc ~ pel p- We can express this in terms of the ef- 
fective mass of the free n e utrons , m*. Then we have (see 
IPrix. Comer fc Anderssoiil (|2002!) for a discussion of the 
analogous problem in a superfluid neutron star core) 



En = 1 



mn 



(35) 
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and it follows that 



rrin 



X 



(36) 



It is easy to show that x'^^"^ encodes the difference between 
the superfluid result and the standard result for a single com- 
ponent crust, i.e. with pc ^ p in (|2ip . This may be the most 
meaningful comparison to make, since all previous studies of 
crust oscillations have assumed the single component model. 

What do we learn from these results? First of all, we 
see that in the limit m* rrin, when the medium effects 
that lead to the effective mass differing from the bare mass 
are not so great, we have — > 1 and x ~* ^c- The waves 
in such a system, cf. (|33p. are the usual shear waves with 
a (as we will see later) relatively small magnetic correction. 
Of course, the results could still differ significantly from the 
standard single component model. The largest effect that 
one would expect would be, for Xc ~ 0.8 cf. Figure [1] a 
frequency increase by about a factor of 2. However, the ef- 
fective mass is expected to be larger than the bare mass so 
let us consider the opposite limit, which may well apply in 
parts of the neutron star crust [see, for example, IChamell 
l|2006l )]. Then we have m* 2> rnn- Using also a^c ^ 1, we see 
that £♦ « a;;, or X ^ 1- In this limit, it is very difflcult for 
the free neutrons to move relative to the crust. The upshot 
of this is that the waves in the system tend to the frequency 
predicted for a pure elastic crust without a superfluid com- 
ponent. 

These two extremes show that the presence of the su- 
perfluid in the neutron star crust can have a significant effect 
on the frequencies of waves in the system. According to the 
data in Figure [TJ the combined effect is at the 10% level 
(compared to the single component crust result). The re- 
sults clearly show that the superfluid component must be 
considered if we want to develop high precision magnetar 
crust seismology. Of course, in reality we are mainly inter- 
ested in the global oscillations. Then the local effects that we 
have worked out will be (in some sense) averaged through- 
out the crust. One may expect this to decrease the role of 
the superfluid since the effective neutron mass may only be 
large in parts of the crust. Of course, the real answer re- 
quires a detailed mode calculation. This problem remains to 
be solved. In order to provide reliable results, such an effort 
should draw on more complete studies of the entrainment 
for the crust superfluid. One should also worry about the 
relevance of vortex pinning and the mutual friction. 

Finally, let us discuss the relative importance of the 
magnetic field. Scaling to "typical" values, we have 



10^ 



108cm/s/ VlOi4g/cm' 



dyne/cm (37) 



Then, it follows from (1331) that we need to consider 



UJA_ 



0.08 



10*cm/s 



B 



10i4g/cmV VlQi^G, 
(38) 

This shows that we can safely ignore the magnetic effects in 
the high-density region of the crust, cf. Figure [1] In order for 
the magnetic term to dominate at the base of the crust we 
need B ~ 10^^ G, stronger than the field strength inferred 
for magnetars. Of course, one has to be a little bit careful 
here. First of all, it is entirely possible that the interior field 



1.5 



10* 



P [g/cm 



10* 




20 



0.001 0.01 
n. [fm ] 



Figure 1. This figure illustrates the density dependence of the 
different parameters that affect the wave propagation in the crust. 
We show, as functions of the total baryon number density , the 
superfluid neutron fraction Xn = rtn/rtb (dash-dot, left scal e) for 
the equation of state discussed bv lDouchin fc Haensel ||2001^ . a fit 
to the effective neutron mass m^/rriu ( solid, right scale) based on 
the numerical results of lChamell | |2006| ) (data points indicated by 
diamonds), and the ratio between the Alfven and the shear wave 
speeds (vj^/vg)^ (dashed, left scale). The overall effect that the 
presence of the crust superfluid has on the local wave propagation, 
compared to the standard single component crust, is (as discussed 
in the main text) represented by x"^''^ (solid, left scale). The 
horizontal dashed grey line indicates unity on the left scale. 



is stronger than the exterior dipole field which leads to the 
observed braking of the magnetar spin. Secondly, (|38p indi- 
cates that the magnetic terms will dominate as we approach 
the surface of the star. However, our analysis breaks down 
completely in the surface region. Basically, the MHD ap- 
proximation is only valid as long as the Alfven wave speed 
is significantly below the speed of light. In order for this to 
be the case, we require 



p> 10** 



B 



1015 G 



(39) 



When this is no longer true, as it will not be in a transi- 
tion region near the surface of any neutron star, one must 
consider the complete Maxwell equations. It is interesting 
we compare (|39|) to the estimated density for t he top of the 
crust, cf. for example Eq. (1) from IPirol (j2005l ). 



ptop ~ 2.3 X 10^ 



3 X 108 K 



where A is the number of nucleons and Z is the charge per 
ion, respectively. This density is scaled to iron, which means 
that the fiducial values should be relevant near the top of 
the crust. This suggests that MHD remains a valid approxi- 
mation throughout a magnetar crust. The transition to elec- 
tromagnetism takes place in the neutron stars envelope. The 
nature of this transition is an important problem that de- 
serves more attention. 
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3 MULTIFLUID DYNAMICS OF THE CORE 

3.1 Lagrangian perturbation equations 

Tlie magnetoliydrodynamics in tlie core is, in tiie simplest 
case, formulated in terms of three distinct fluids associated 
with the neutrons, protons and electrons. The former two 
particle species are expected to be superfluid and supercon- 
ducting, respectively. Due to the smallness of the electron 
mass, the electron fluid degree of freedom can be suppressed 
and the MHD equations effectively lead to a two-fluid model 
l|Mendelllll993) . The protons in the outer core arc expected 
to fo rm a type II superconductor (Bavm, Pethi ck & Pinesj 
which means that the magnetic fleld is carried by 
a large number of fluxtubes, each with a flux quantum 
(po ~ hc/2e. This should be the case provided the mag- 



netic fleld is below a critical value Hc2 



3 X 10^" G 



l|Bavm. Pethick fc PinesI Il969l ). This critical value repre- 
sents the fleld strength at which the magnetic fluxtubes over- 
lap and can no longer be treated individually. Above this 
threshold the magnetic field behaves "classically". It should 
be noted that, even though the critical fleld is large, this is 
not an unrealistic possibility for magnetars given that the 
magnetic field in the interior could be considerably higher 
than the exterior dipole fleld. Hence, one ought to consider 
both superconducting and normal protons. The latter case 
is, however, trivial. The desired result follows immediately 
from the previous section, e.g. (|33p . if we take the limit 

For a non-rotating star, we again neglect the vortex- 
mediated mutual friction and the neutron vortex ten- 
sion. Omitting also a small entrainment induced mag- 
netic term that originates from the London field by 
means of which the proton supercon d uctor rotates [see 
iGlampedaJtis, Samuelsson fc AnderssonI l|2008l ) for discus- 
sion], the superfluid neutron dynamics is still governed by 
©• 

The combined proton-electron dynamics is 
a little bit more complicated. As dis cussed by 
iGlampedakis, Samuelsson fc AnderssonI l|2008l ') the rel- 
evant equation of motion takes the form 

{dt + viV,) K + ecwD -f V,(/i + $) + e.w^'V.vi 

^-(f^ + t^) (41) 

where vl and are the velocites and ~ Vn — 

vl, and we have (again) neglected the mutual friction. 
There are two force terms on the right-hand side of 
this equation. The first, fl, is the usual electromagnetic 
Lorentz force given by ((Sj. The second, tl, represents 
the smooth-averaged tension of the magnetic fluxtubes. 
Remarkably, the Lorentz force does not play a role in 
th e final superconducting MHD equatioii s . As discussed 
by [G lamped akis. Samuelsson fc AnderssonI (|2008l ) [see also 
l|Me ndcll 199^)] it is exactly cancelled by a term originating 
from the fluxtube tension. In the case of a non-rotating star 
this leads to the magnetic force taking the form 



47rc 



[V,(//c 



(42) 



where the lower critical magnetic fleld Hci — Hci{pc) ~ 
10^^ G. 

We have deflned the vector Wc, representing the (av- 



eraged) canonical proton vorticity l|Prixll2005l ). This means 
that we havcQ 

Wl ^ —B' (43) 
nicC 

For a uniform incompressible model, this means that 

In order to close the system of the MHD equations we 
need to provide a relation between the magnetic field and 
the fluid velocity. This relation follows from the mag- 
netic induction equation. Neglecting the coupling forces be- 
tween the electrons and the neutron an d proto n fluids (see 
iGlampedakis. Samuelsson fc AnderssonI l|2008l ) for discus- 
sion) the induction equation takes the standard form. 



(45) 

Its perturbed form is given by (|10|) from which we obtain 
the Eulerian perturbation of the magnetic fleld. 



(46) 



Using this result, we flnd (for a uniform background and 
incompressible perturbations) 

^ cl&B'ig'" - B'B'^)V,V,Cfe (47) 



where 



2 HciB 

ca = 

47rpc 



(48) 



We now have all the relations we need to discuss short wave- 
length waves in the superfiuid/superconducting system. 

3.2 Plane- wave analysis 

Most of the analysis works out exactly as in the crust prob- 
lem (obviously in the fi ^ limit). The only difference is the 
form of the perturbed magnetic force. With the plane-wave 
assumption, we see that 



(49) 



Comparing the magnetic forces in the normal and supercon- 
ducting cases (equations (fM)l and (fl5)) . respectively) we note 
some important differences. The characteristic speeds are 
obviou sly different, cj = (Hd/ B)v\. This is a well-known 
effect (jEasson fc Petliicklll977l : lMendelllll998l) . The two ve- 
locities would differ by a factor ~ 10'^ for a canonical pulsar 
with B = 10^^ G. However, the effect will not be as dramatic 
for magnetars for which B ^ Hd « 10^^ G. Another differ- 
ence arises from the second term in (|49p . This term is present 
provided B'-Q / 0, or in other words when the wave- vector 
fc' does not point along the direction of the magnetic field. 



^ Strictly speaking, this expression is only valid for the non- 
rotating background. Perturbations of will also contain the 
perturbed fluid velocities. However, as long as we are focussing 
on the leading order contributions these can be neglected. That 
this is a legitimate approximation is easy to see since the charac- 
teristic frequency 

B Si 10^^ — — 

■mcc VlO^^G^ 

is much higher than any other relevant frequency in the problem. 
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As a result, the magnetic force for superconducting protons 
has a non-zero component along the magnetic field. Mean- 
while, we have _Bi Ac/mag = in the normal proton case, cf. 
H24p . In the plane- wave problem the terms that give rise to 
this difference vanish identically. However, this rather subtle 
difference in the magnetic forces will be relevant in many less 
idealised situations, e.g. for global mode oscillations. Hence, 
one must be careful before using intuition gained from stan- 
dard MHD problems in the case of a superconducting core. 
There is certainly mor^fl to the problem than a simple "re- 
placement" BHci . 

Combining the perturbation equations as in the previ- 
ous section, we readily arrive at 



Cj^k (kjlB- 



Al = ~clk'{k,B'f{B,Ai)B^ (50) 



Working things out as in the crust case, we project this 
equation onto ki. Since kiAl — it then follows that we 
have the constraint 



ciCkjB'fiBjAi)=0 



(51) 



This has the same implications as in the crust problem. It 
immediately follows that the right-hand side of (|5Up must 
vanish. Hence, we have the dispersion relation. 



£i,Cj^k {kjB 



47rpc 



(52) 



Our main interest here concerns the role of the super- 
fluid neutron component. Its presence is reflected by the 
entrainment factor in (|52[) . To quantify its relevance, we ex- 
press the entrainment in terms of the effective proton mass, 
i.e. we use £c = l — rrip/m-p. Then it follows, since the proton 
fraction in the core is small, that 



^dk' 



(53) 



Since it is expected that 0.3 < m'/mp < 0.7 (see 
IPrix. Comer fc AnderssonI (|2002l ) for discussion) we see that 
the presence of the superfluid neutrons will lead to a ~ 
20 — 80% increase in the frequency of the core waves. This 
effect is large enough that cannot be neglected. It may, in 
fact, also be observable. If one accepts the argument that the 
magnetic field couples motion in the crust to the core, and 
that the core fluid is therefore partaking on the oscillation, 
then the entrainment will affect the observed frequencies. 

We cannot at this point say much about the global os- 
cillations of a magnetic neutron star core; it is a problem 
that remains to be solved in detail. It is co mplicated by the 
likely presence of an "Alfven continuum" (|Levinl 12007 ). At 
this point it is not clear to what extent the continuum pre- 
vails in more detailed neutron star models. However, it is 
easy to see how the presence of the superfluid will mani- 
fest it self in the continuum toy-model considered bv iLeviiil 
l|2007l ). The frequency range of the continuum will simply 
scale according to (f52]). 



^ For an interesting recent discussion on how superconduc- 
tivity may affect the sta bility properties of the star, see 



lAkgun fc WassermanI | |200^ . 



4 CONCLUDING REMARKS 

In this paper we have investigated the role of neutron 
star superfluidity for magnetar oscillations. The results im- 
pact on attempts to use data from observed quasiperi- 
odic oscillations in the tails of magnetar flares to place 
contraints on neutron s t ar pa rameters, see for example 
ISamuelsson fc AnderssonI (|2007l ). Using a plane- wave anal- 
ysis we estimated the effects of the neutron superfluid in 
the elastic crust region. This, the flrst ever, analysis of the 
combined magnetic-elastic-superfluid crust problem demon- 
strated that the superfluid imprint is likely to be more sig- 
nificant than the effects of the crustal magnetic field. This 
is, of course, assuming that the SGR flare mechanism does 
not deposit sufficient heat in the crust to raise the system 
above the superfluid transition tem perature. Available esti- 
mates, e.g. iKouveliotou et all l|2003h . suggest that this is un- 
likely. We also considered the region immediately beneath 
the crust, where superfluid neutrons are thought to coexist 
with a type II proton superconductor. Since the magnetic 
field in the latter is carried by an array of fl uxtubes, see 
IClampedakis. Samuelsson fc AnderssonI (|2008l ) for a recent 
discussion, the dynamics of this region differs from standard 
magnetohydrodynamics. We showed that the presence of the 
neutron superfluid (again) affects the oscillations of the sys- 
tem. This acco rds well with previous results of, in particular, 
iMendelll l|l998l ). 

Our estimates show that the superfluid components 
cannot be ignored in efforts to carry out magnetar seismol- 
ogy. This increases the level of complexity of the modelling 
problem, but also points to the exciting possibility of using 
observations to probe the superfluid nature of supranuclear 
matter. Future work needs to extend our analysis to consider 
the global oscillations of magnetic-superfluid-elastic neutron 
stars. This is a very interesting problem because, in addi- 
tion to enabling a more detailed seimsology analysis, it may 
also provide insight into ro tational glitches in magnetars 
(jPib. Kaspi fc Gavrii]|[l993 ). It is generally believed that 
superfluidity plays a key role in radio pulsar glitches. Re- 
cent dev elopments in modelling these even ts is showing some 
promise (|Glampedakis fc Andersso"nll2008l ) , and it would ob- 
viously be highly relevant to extend this analysis to strongly 
magnetised systems. 
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